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We describe searches for B meson decays to the charmless vector- vector final states uK* , up, uuj, 
and lj4> with 233x10® BB pairs produced in e"'"e~ annihilation at y/s = 10.58 GeV. We also search 
for the vector-scalar B decay to ujfo. We measure the following branching fractions in units of 10~^: 
^ ioK*'') = 2.4 ± 1.1 ± 0.7 (< 4.2), B{B+ ^ luK*+) = 0.6t J] (< 3.4), Z3(B'' ^ V) = 
-0.6 ± O.ltol (< 1.5), B{B+ Lup+) = 10.6 ± 2.lli;f,, B{B° luuj) = l.SlJJ ± 0.4 (< 4.0), 
B{B° LO(l>) = 0.1 ± 0.5 ± 0.1 (< 1.2), and B{B° -> ujfo) = 0.9 ± 0.41°;^ (< 1.5). In each case 
the first error quoted is statistical, the second systematic, and the upper limits are defined at the 
90% confidence level. For B"*" —> tjp"*" decays we also measure the longitudinal spin component 
/l = 0.82 ± 0.11 ± 0.02 and the charge asymmetry Acp = 0.04 ± 0.18 ± 0.02. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Until recently, hadronic decays of B mesons to pairs 
of light vector mesons (VV final states) have received 
less theoretical and experimental attention than decays 
to two pseudoscalar mesons (PP) or one pseudoscalar 
and one vector meson (PV). Early papers presented cal- 
culations for branching fractions, CP-violating asymme- 
tries , and relative spin component contributions for 
these decays. The measurement three years ago of an un- 
expectedly small value of the fraction of the longitudinal 
spin component (/i) in penguin-dominated B —f (pK* 
decays 1^ triggered new theoretical activity. There 
have been several attempts to understand the small value 
of Jl within the Standard Model (SM)j5| and many pa- 
pers suggested non-SM explanations |3- Further infor- 
mation about these effects can come from both branch- 
ing fraction and measurements in decays such as 
B ujK* or B uj(p, which are conjugate to B ^ (f)K* 
via an SU(3) rotation Q. Information on these and re- 
lated charmless decays can additionally be used to pro- 
vide sensitivity to the CKM angles a and 7 0. 



We have discussed above mostly penguin-dominated 
decays. There are also decays with the K* replaced by 
a p, w or meson where tree diagrams are expected to 
be more important. These include B decays to the final 
states PP, up, and louj. The decay B on is known to 
be nearly fully longitudinally polarized [9|, [lOj and these 
other predominantly tree decays are expected to behave 
similarly j with branching fractions as predicted in Q . 



We report results of measurements of B decays to the 
charmless VV final states loV , where V represents a neu- 
tral or charged K* or p, or an a; or meson. We also 
measure the decay B^ — > a;/o(980) which shares the same 
final state as the B^ ujp^ decay. Due to the current 
small signal samples, only the branching fractions and 
the fraction of the longitudinal spin component are mea- 
sured, the latter by integrating over the azimuthal angles, 
for which the azimuthal acceptance is uniform. The an- 



gular distribution is 



r c? cos 01 d cos 62 

I 1 1 (1 - /l) sin^ 01 sin2 02 + /l cos^ 0i cos^ 62 1 , 

where 0^ is the helicity angle in the 14 rest frame 
with respect to the boost axis from the B rest frame 
and /i is the fraction of the longitudinal spin com- 
ponent. For B^ —> ujp~^, we also measure the direct 
CP-violating time-integrated charge asymmetry Acp = 
(r^ — r"'")/(r~ -l- r"*"), where the superscript on the F 
corresponds to the sign of the B^ meson. 

The results presented here are based on data collected 
with the BABAR detector I?] at the PEP-II asymmet- 
ric e+e" collider located at the Stanford Linear Accel- 
erator Center. An integrated luminosity of 211 fb~^, 
corresponding to 232.8x10^ BB pairs, was recorded 
at the T(4S') resonance (center-of-mass energy ^/s = 
10.58 GeV). 

Charged particles from the e^e~ interactions are de- 
tected, and their momenta measured, by five layers of 
double-sided silicon microstrip detectors surrounded by 
a 40-layer drift chamber, both operating in the 1.5-T 
magnetic field of a superconducting solenoid. We iden- 
tify photons and electrons using a CsI(Tl) electromag- 
netic calorimeter (EMC). Further charged particle iden- 
tification (PID) is provided by the average energy loss 
(dE / dx) in the tracking devices and by an internally re- 
flecting ring-imaging Cherenkov detector (DIRC) cover- 
ing the central region. 

We reconstruct the P-daughter candidates through 
their decays p° tt+tt^, /o(980) tt+tt^, p+ 
tt+tt", X+7r-(denoted K*j^+ _), K*+ 

K+ir^{K*+^^o), K*+ ^ X07r+(if^t^J,"c.^^+7r-^o, 
(f) K^K", 7r° 77, and K^^ tt+tt^ (charge- 
conjugate decay modes are implied throughout). Table 
^ lists the requirements on the invariant mass of these 
particles' final states. For the p, K* , (j), and lo invari- 
ant masses these requirements are set loose enough to 
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TABLE L Selection requirements on the invariant mass and 
helicity angle of B-daughter resonances. The helicity angle is 
unrestricted unless indicated otherwise. 



State 



inv. mass ( MeV) 



helicity angle 



K'+o^+ 755 <mK^< 1035 -0.85 < cos e < 1.0 
^k+ttO 755 <mK^< 1035 -0.8 < cos 61 < 1.0 

0.85 < cose < 0.85 



pVfo 
P+ 



755 < niKn < 1035 
540 < m^^ < 1060 
470 < m^^ < 1070 
740 < m-jr-KTT < 820 
1009 < rriKK < 1029 
120 < m^j < 150 
473 < m^^ < 522 



-0.7 < COS0 < 0.85 



include sidebands, as these mass values are treated as ob- 
servables in the maximum-likelihood fit described below. 
The p° and /o (we use /o as shorthand for /o(980)) yields 
are extracted from a simultaneous fit to the same data 
sample. For candidates we further require the three- 
dimensional flight distance from the event primary vertex 
to be greater than three times its uncertainty. Secondary 
pion and kaon candidates in p, K* , and uj candidates are 
rejected if their DIRC, dE/dx, and EMC PID signature 
satisfies tight consistency with protons or electrons, and 
the kaons (pious) must (must not) have a kaon signature. 

Tabled also gives the restrictions on the K* and p he- 
licity angles 9 (previously defined for Eq. ^ , imposed to 
avoid regions of rapid acceptance variation or combina- 
torial background from soft particles. To calculate 9 we 
take the angle relative to, for uj, the helicity axis of the 
normal to the decay plane, for p and 0, the positively- 
charged (or only charged) daughter momentum, and for 
K* the daughter kaon momentum. 

A i3-meson candidate is characterized kinemat- 
ically by the energy-substituted 

and 



mass 



the 



energy differ- 



ence AE = E^ — \-/s, where (£'o,po) and {Eb,Pb) 
are four-momenta of the T(45') and the B candidate, 
respectively, s is the square of the center of mass 
energy, and the asterisk denotes the T{4S) frame. The 
resolution on A^; (toes) is about 30 MeV (3.0 MeV). 
Our signal falls in the region |Ai?| < 0.1 GeV and 
5.27 < Toes < 5.29 GeV, which is then extended to 
include sidebands that provide a good description of the 
background. The average number of candidates found 
per selected event in data is in the range 1.1 to 1.7, 
depending on the final state. We choose the candidate 
with the smallest value of a constructed from the 
deviations of the B daughter resonance masses (all 
particles in Table HI except 7r°, if", and /o) from their 
expected |0 values. 

Backgrounds arise primarily from random combina- 
tions of particles in continuum e~^e~ — > qq events {q = 
u,d,s,c). We reduce these by using the angle 9t be- 



tween the thrust axis of the B candidate in the T{4S) 
frame and that of the rest of the charged and neutral 
particles in the event. The distribution of | cos 0x1 is 
sharply peaked near 1.0 for qq ]et pairs, and nearly uni- 
form for B meson decays. The requirements, chosen to 
reduce the sample size for the large background modes, 
are |cos0t| < 0.9 for the cucj) niode, |cos0t| < 0.8 for 
the ojuj and ujK* modes, and |cos6't| < 0.7 for the ujp 
modes. In the maximum-likelihood fit described below, 
we also use a Fisher discriminant that combines four 
variables defined in the T(4S') frame: the polar angles 
with respect to the beam axis of the B momentum and 
B thrust axis, and the zeroth and second angular mo- 
ments, Lq and L2, of the energy flow about the B thrust 
axis in the T(4S') frame. The moments are defined by 
Lj = J2iPi ^ |cos0i|"' , where 9i is the angle with respect 
to the B thrust axis of a charged or neutral particle i, pi 
is its momentum, and the sum excludes the B candidate's 
daughters. 

From Monte Carlo (MC) simulation 0| we determine 
the most important charmless BB backgrounds (typi- 
cally about a dozen background modes for each signal 
final state). We include a variable yield for these in the 
fit described below. We also introduce a component for 
non-resonant tttt and Kn background. The magnitude 
of this component is fixed in the fit as determined from 
extrapolations from higher-mass regions. 

We obtain yields and values of /l and Acpbom ex- 
tended unbinned maximum-likelihood fits with input ob- 
servables AE, TOes, -?"and, for the vector meson fc, the 
mass TOfe and helicity angle 9k- For each event i and hy- 
pothesis j (signal, qq background, BB background) we 
define the probability density function (PDF) 

V} - V,{mEslV,{AE^)V,{nV,{rnlml,9l9l). (2) 

We check for correlations in the background observables 
beyond those accounted for in this PDF and find them 
to be small. For the signal component, we correct for the 
effect of small neglected correlations (see below). The 
likelihood function is 



Nl 



(3) 



i=l j 



where Yj is the yield of events of hypothesis j and N is 
the number of events in the sample. 

The PDF factor for the resonances in the signal takes 
the form Vi,sig{m\)r2,sig{m'2)Qi9l,9l) with Q given by 
Eq.nmodified to account for detector acceptance. For qq 
background it is given for each resonance independently 
by V,^{ml,9l) = Vp^iml)Vp^i9l) +V.{ml)VM), dis- 
tinguishing between true resonance (T^pk) and combina- 
torial (Vc) components. For the BB background we as- 
sume that all four mass and helicity angle observables are 
independent. 



6 



TABLE II: Signal yield Y and its statistical uncertainty, bias Yq, detection efficiency e, daughter branching fraction product 
H'Bi, significance 5* (with systematic uncertainties included), measured branching fraction B, 90% C.L. upper limit, measured 
or assumed longitudinal polarization, and Acp- 



Mode 


Y Yo 
(events) (events) 


e 
(%) 


(%) 


S 

{^) 


B 

(10-^) 


B U.L. 

(io-«) 




Acp 






11 


13.2 


59.2 


2.4 


2.4 ± 1.1 ±0.7 


4.2 


^- ' -"--0.24 


— 


ujK''+ 


-3.6l^° 


-5 
6 


12.5 
8.0 


20.3 
29.6 


0.1 

0.5 
0.4 


n 9+1.7 +1.5 
^■'^-l.b -1.1 
1 1 -t-^.o-f- l.O 
^•^-2.1-1.2 
n c + 1. 4+1.1 
^^•"-1.2-0.9 


4.2 

5.7 
3.4 


0.7 fixed 

0.7 fixed 
0.7 fixed 


— 






-2 


11.6 


89.1 


0.6 


-0.6 ± 0.71;; j 


1.5 


0.9 fixed 




^fo 




4 


15.2 


59.4 


2.8 


0.9 ± 0.4t°;^ 


1.5 






LOp^ 


1561^? 


11 


6.6 


89.1 


5.7 


10.6±2.lli;^ 




0.82 ±0.11 ±0.02 


0.04 ±0.18 ±0.02 


LOLO 




8 


12.9 


77.5 


2.1 


1-8^0.9 ±0.4 


4.0 


0.79±0.34 




LO(j) 


q 1 +4.4 
O-J^-8.5 


1 


19.0 


43.2 


0.3 


0.1 ±0.5 ±0.1 


1.2 


0.88 fixed 





For the signal, BB background, and non-resonant 
background components we determine the PDF param- 
eters from simulation. We study large control samples 
of i? — > Dt: decays of similar topology to verify the sim- 
ulated resolutions in A£' and mEs, adjusting the PDFs 
to account for any differences found. For the continuum 
background we use (toes ^ Ai?) sideband data to obtain 
initial values, before applying the fit to data in the signal 
region, and leave them free to vary in the final fit. 

The parameters that are allowed to vary in the fit in- 
clude the signal, BB background, and non-resonant back- 
ground yields, and continuum background PDF param- 
eters. For the three modes with signals of more than 2cr 
significance, we vary in the fit to properly account for 
the variation of efficiency with f^. For B^ ujp'^ we 
also vary the signal and background charge asymmetries. 
For the fits with little signal, we fix to a value that is 
consistent with a priori expectations (see Table llT|) . and 
account for the associated systematic uncertainty. 

To describe the PDFs, we use the sum of two Gaussians 
for 7^sig(TOEs)5 'Psigi'^E), and the true resonance compo- 
nents of Vj{mk)', for we use an asymmetric Gaussian 
for signal with a small Gaussian component for back- 
ground to account for important tails in the signal re- 
gion. The background toes shape is described by the 
function 1 — exp [—^(1 — x^)] (with x = toes/^'b) 
and the distributions of masses nik by second or third or- 
der polynomials. The background PDF parameters that 
are allowed to vary in the fit are ^ for toes, slope for AE, 
the polynomial describing the combinatorial component 
for TOfc, and the peak position and lower and upper width 
parameters for J^. 

We evaluate possible biases from our neglect of corre- 
lations among discriminating variables in the PDFs by 
fitting ensembles of simulated experiments into which 
we have embedded the expected number of signal events 
and BB background events, randomly extracted from the 
fully-simulated MC samples. We give in Table the 



yield bias Iq found for each mode. Since events from 
a weighted mixture of simulated BB background decays 
are included, the bias we measure includes the effect of 
this background. 

The systematic uncertainties on the branching frac- 
tions arising from lack of knowledge of the PDFs have 
been included in part in the statistical error since most 
background parameters are free in the fit. For the sig- 
nal, the uncertainties in PDF parameters are estimated 
from the consistency of fits to MC and data in large con- 
trol samples of topology similar to signal. Varying the 
signal PDF parameters within these errors, we estimate 
yield uncertainties for each mode. The uncertainty in the 
yield bias correction is taken to be the quadratic sum of 
two terms: half the bias correction and the statistical 
uncertainty on the bias itself. Similarly, we estimate the 
uncertainty from modeling of the BB backgrounds as the 
change in the signal yield when the number of fitted BB 
events is fixed to be within one sigma of the expected 
number of BB events from MC. For the non-resonant 
TTTT or Kt: backgrounds the uncertainty is taken as the 
change in the signal when the background yield is var- 
ied within the uncertainty of the fits to the higher-mass 
regions. For modes with fixed /l, the uncertainty due 
to the dependence of signal efficiency on /l is evaluated 
as the measured change in branching fraction when is 
varied by ±0.3 (up to a maximum of /l = 1). These ad- 
ditive systematic errors are dominant for all modes; the 
PDF variation is always the smallest but the others are 
typically similar in size. 

Uncertainties in our knowledge of the efficiency, found 
from studies of data control samples, include 0.8% x Nt, 
3.0% X iV^o , and 1% for a decay, where Nt is the num- 
ber of tracks, and N^o is the number of 7r°s in a decay. 
We estimate the uncertainty in the number of B mesons 
to be 1.1%. Published data provide the uncertainties 
in the i3-daughter product branching fractions (1-2%). 
The uncertainties from the event selection are 1-2% for 
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the requirement on cos 9t . 

The systematic uncertainty on for —> uip^ in- 
cludes the effects of fit bias, PDF-parameter variation, 
and BE and non-resonant backgrounds, all estimated 
with the same method as for the yield uncertainties de- 
scribed above. From large inclusive kaon and B-decay 
samples, we find a systematic uncertainty for Acp of 0.02 
due mainly to the dependence of reconstruction efficiency 
on the momentum of the charged p daughter. We find 
for the B+ ujp+ background, A'^Jp = -0.010 ± 0.007, 
confirming this estimate. 

In Table ITU we also show for each decay mode the mea- 
sured branching fraction with its uncertainty and signif- 
icance together with the quantities entering into these 
computations. The significance is taken as the square 
root of the difference between the value of —2 In £ (with 
systematic uncertainties included) for zero signal and 
the value at its minimum. For all modes except for 
5+ — > ujp^ we quote a 90% confidence level (C.L.) up- 
per limit, taken to be the branching fraction below which 
lies 90% of the total of the likelihood integral in the pos- 
itive branching fraction region. In calculating branching 
fractions we assume that the decay rates of the T{AS) to 
B+B- and B°B° are equal For decays with K*+, 

we combine the results from the two K* decay channels, 
by adding their values of —2 In £, taking into account the 
correlated and uncorrelated systematic errors. 

We present in Fig. the data and PDFs projected 
onto mES and AE, for subsamples enriched with a mode- 
dependent threshold requirement on the signal likelihood 
(computed without the PDF associated with the variable 
plotted) chosen to optimize the significance of signal in 
the resulting subsample. Fig. |21 gives projections of cos 6 
for the B~^ — > ojp'^ decay. 

The branching fraction value B given in Table ^ for 
_B+ Lup'^ comes from a direct fit with the free param- 
eters B and /l, as well as Acp- This choice exploits the 
feature that B is less correlated with than is either 
the yield or efficiency taken separately. The behavior of 
-2ln C{fL,B) is shown in Fig.Ol 

In summary, we have searched for seven charmless 
hadronic _B-meson decays. We observe -B^ Lup'^ with a 
significance of 5.7 a, and establish improved 90% C.L. up- 
per limits for the other modes, with the following branch- 
ing fractions: 



B{B°^ 




= 2.4 ± 1.1 ±0.7 (< 4.2) 


X 10" 


6 


B{B+ 


ujK*+) 


= O.ellit^i (< 3.4) X 10-^ 




BiB"" 




= -0.6±0.7^°j (< 1.5) 


X 10" 


6 


B{B+ - 


-.UP+) 


= 10.6±2.lt} J^ X 10"^ 






B{B° 


LOLO) 


= 1.8to J ±0.4 (< 4-0) X 


10"^ 




B{B" 


UJ(j)) 


= 0.1 ±0.5 ±0.1 (< 1.2) 


X 10" 


6 


B{B° 




= 0.9 ±0.41^-2 (< 1.5) X 


10"^ 





In each case the first error quoted is statistical, the second 




-0.2 -0.1 0.1 0.2 5.25 5.26 5.27 5.28 5.29 



A E (GeV) Mjjs (GeV/c^) 

FIG. 1: Projections of AE (left) and tties (right) of events 
passing a sig nal likelihood threshold for B° ^ uK*° , (a,b), 
B+ ujK*+, (c,d) B° up°, (e,f) B+ ujp+, (g,h) 
B° ujuj, (i,j) B° ^ uj^, (k,l) and B° ufo, (m,n). The 
solid curve is the fit function, the dashed curve is the sig- 
nal contribution, and the dot-dashed curve is the background 
contribution. 




IcosO"! coseP 



FIG. 2: Projections of helicity- angle cosines, of events pass- 
ing a signal likelihood threshold for to (left) and p"*" (right) 
from the fit for i?"*" — > a;p^ decay. The solid curve is the fit 
function, the dashed curve is the signal contribution, and the 
dot-dashed curve is the background contribution. 
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Ql 1 1 1 1 1 

5 10 15 20 25 

Branching Fraction(B"^^ CO p"^)(xlO'*) 

FIG. 3: Distribution of -2ln £{fL,B) for B+ ujp+ decay. 
The solid dot gives the central value; curves give the contours 
in 1-sigma steps (A^— 2 In S) = 1). 

systematic, and the upper limits are taken at 90% C.L. 

For B'^ — > Lvp'^ we also measure the longitudinal 
spin component /l = 0.82 ± 0.11 ± 0.02 and charge 
asymmetry Acp = 0.04 ± 0.18 ± 0.02. The longitudi- 
nal spin component is dominant, as it is for B ^ pp 
|9|, [lOj . Assuming tree dominance we would naively ex- 
pect the branching fraction for B^ ojp'^ to be equal 
to that of B^ P^P^- However the measured branch- 
ing fraction for B{B^ '-^P'^) is more than two stan- 
dard deviations smaller than the world average B{B^ 
p+p°)=26±6 X 10-s [13. 

Our branching fraction results are in agreement with 
theoretical estimates [HllH. 
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